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Abstract — In this work, we introduce convolutional codes for 
network-error correction in the context of coherent network 
coding. We give a construction of convolutional codes that correct 
a given set of error patterns, as long as consecutive errors are 
separated by a certain interval. We also give some bounds on the 
field size and the number of errors that can get corrected in a 
certain interval. Compared to previous network error correction 
schemes, using convolutional codes is seen to have advantages in 
field size and decoding technique. Some examples are discussed 
which illustrate the several possible situations that arise in this 
context. 

I. Introduction 

Network coding was introduced in [1] as a means to 
improve the rate of transmission in networks, and often achieve 
capacity in the case of single source networks. Linear network 
coding was introduced in [2]. network-error correction, which 
involved a trade-off between the rate of transmission and the 
number of correctable network-edge errors, was introduced in 
[5] as an extension of classical error correction to a general 
network setting. Along with subsequent works [6] and [7], 
this generalized the classical notions of the Hamming weight, 
Hamming distance, minimum distance and various classical 
error control coding bounds to their network counterparts. An 
algebraic formulation of network coding was discussed in [3] 
for both instantaneous networks and networks with delays. In 
all of these works, it is assumed that the sinks and the source 
know the network topology and the network code, which is 
referred to as coherent network coding. 

Random network coding, introduced in [4] presented a 
distributed network coding scheme where nodes independently 
chose random coefficients (from a finite field) for the linear 
mixing of their inputs. Subspace codes and rank metric codes 
were constructed for the setting of random network coding in 
[8] and [9]. 

Convolutional network codes were discussed in [10]-[12] 
and a connection between network coding and convolutional 
coding was analyzed in [13]. In this work, convolutional 
coding is introduced to achieve network-error correction. We 
assume an acyclic, single source, instantaneous (delay-free) 
network with coherent linear network coding to multicast 
information to several sinks. 

We define a network use as a single usage of all the edges 
of the network to multicast utmost min-cut number of symbols 
to each of the sinks. An error pattern is a subset of the set of 
edges of the network which are in error. It is seen that when the 
source implements a convolutional code to send information 
into the network, every sink sees a different convolutional 
code. We address the following problem. 



Given an acyclic, delay-free, single-source network with a 
linear multicast network code, and a set of error patterns $, 
how to design a convolutional code at the source which shall 
correct the errors corresponding to the error patterns in $, as 
long as consecutive errors are separated by a certain number 
of network uses? 

The main contributions of this paper are as follows. 

• For networks with a specified network code, convolu- 
tional codes have been employed to achieve network-error 
correction for the first time. 

• An explicit convolutional code construction (for the net- 
work with a given network code) that corrects a given 
pattern of network-errors (provided that the occurrence 
of consecutive errors are separated by certain number of 
network uses) is given. 

• The convolutional codes constructed in this paper are 
found to offer certain advantages in field size and de- 
coding over the previous approaches of block network- 
error correction codes (BNECCs) of [7] for network error 
correction. 

• Some bounds are derived on the minimum field size 
required, and on the minimum number of network uses 
that two error events must be separated by in order that 
they get corrected. 

The rest of the paper is organized as follows. Section |ll] 
gives a primer on convolutional codes and MDS convolu- 
tional codes. In Section Hill we discuss the general network 
coding set-up and network-errors. In Section IIVI we give a 
construction for a input convolutional code which shall correct 
errors corresponding to a given set of error patterns. In Section 
[Vl we give some examples for this construction. In Section 
IVII we compare the advantages and disadvantages of our 
network error correcting codes with that of [7]. Finally, a 
short discussion on the construction of Section HV] constitutes 
Section IVIII along with several directions for further research. 



II. Convolutional codes-Basic Results 

In this section, we review the basic concepts related to 
convolutional codes, used extensively throughout the rest of 
the paper. For g, power of a prime, let Fg denote the finite 
field with q elements. 

For a convolutional code, the information sequence u = 
[uo,Mi, ...,Mt] [ui G F^) and the codeword sequence (output 
sequence) v = [vq, f i, ft] {vi e F^) can be represented in 



terms of the delay parameter z as 

t t 
u{z) ~ UiZ^ and v(z) ~ Viz^ 

i=0 1=0 

Definition 1 ( [14]): A convolutional code, C of rate 
b/c {b < c) is defined as 

C={v{z)e¥^^[[z]] I v{z) ^ u{z)Giz)} 

where G{z) is a 6 x c generator matrix with entries from 
Fg(2)(the field of rationals functions over Fg) and rank b over 
Fg(z), and v{z) being the code sequence arising from the 
information sequence, u{z) e F^[[z]], the set of all 6- tuples 
with elements from the formal power series ring F(j[[z]] over 
F,. 

Two generator matrices are said to be equivalent if they 
encode the same convolutional code. A polynomial generator 
matrix [14] for a convolutional code C is a generator matrix 
for C with all its entries from Fg[z], the ring of polynomials 
over ¥q. It is known that every convolutional code has a 
polynomial generator matrix [14]. Also, a generator matrix 
for a convolutional code is catastrophic [ 14] if there exists an 
information sequence with infinitely many non-zero compo- 
nents, that results in a codeword with only finitely many non- 
zero components. For a polynomial generator matrix G{z), 
let gij{z) be the element of G{z) in the i*^ row and the j'*'* 
column, and Vi :— maxj deg{gij{z)) be the i*^ row degree of 
G{z). Let 5 := J2^i=i be the degree of G{z). 

Definition 2 ( [14] ): A polynomial generator matrix is 
called basic if it has a polynomial right inverse. It is called 
minimal if its degree 5 is minimum among all generator 
matrices of C. 

Forney in [15] showed that the ordered set {j/i, 1^2, I'h} 
of row degrees (indices) is the same for all minimal basic 
generator matrices of C (which are all equivalent to one 
another). Therefore the ordered row degrees and the degree 
5 can be defined for a convolutional code C. A rate h/c 
convolutional code with degree 6 will henceforth be referred 
to as a (c, 6, 5) code. Also, any minimal basic generator matrix 
for a convolutional code is non-catastrophic. 

Definition 3 ( [14] ): A convolutional encoder is a physical 
realization of a generator matrix by a linear sequential circuit. 
Two encoders are said to be equivalent encoders if they encode 
the same code. A minimal encoder is an encoder with minimal 
delay elements among all equivalent encoders. 

Definition 4 ( [14]): The free distance of the convolutional 
code C is given as 

dfree{C) = min {w H {v {z))\v (z) e C,v{z) ^ 0} 

where wh indicates the Hamming weight over F^ . 

A. MDS convolutional codes 

In this subsection, we discuss some results on the existence 
and construction of Maximum Distance Separable (MDS) 
convolutional codes. In Subsection IIV-EI we use these results 
to obtain some bounds on the field size and the error correcting 



capabilities of such MDS convolutional codes when they are 
used for network-error correction. The following bound on the 
free distance, and the existence of codes meeting the bound, 
called MDS convolutional codes, was proved in [16]. 

Theorem 1 ( [16]): For every base field F and every rate 
k/n convolutional code C of degree 5, the free distance is 
bounded as 

dfreeiC) < {fl ~ k){[S / k\ + 1) + S + 1. 

Theorem [T] is known as the generalized Singleton bound. 

Theorem 2 ( [16]): For any positive integers k < n, S and 
for any prime p there exists a field Fg of characteristic p, and 
a rate k/n convolutional code C of degree 6 over Fg, whose 
free distance meets the generalized Singleton bound. 

A method of constructing MDS convolutional codes based 
on the connection between quasi-cycUc codes and convolu- 
tional codes was given in [17]. It is known [17] that the field 
size q required for a (n, k, S) MDS convolutional code C in 
the construction in [17] should be a prime power such that 
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n\{q-l) eindq>d-^^——^+2. (1) 

III. Convolutional Codes for network-error 
Correction - Problem Formulation 

A. Network model 

We consider only acyclic networks in this paper the model 
for which is as in [11]. An acyclic network can be represented 
as a acyclic directed multi-graph Q = (V,£) where V is the 
set of all vertices and £ is the set of all edges in the network. 

We assume that every edge in the directed multi-graph 
representing the network has unit capacity (can carry utmost 
one symbol from Fg). Network links with capacities greater 
than unit are modeled as parallel edges. The network is 
assumed to be instantaneous, i.e, all nodes process the same 
generation (the set of symbols generated at the source at a 
particular time instant) of input symbols to the network in a 
given coding order (ancestral order [11]). 

Let s e V be the source node and T be the set of all 
receivers. Let be the unicast capacity for a sink node T ^ T 
i.e the maximum number of edge-disjoint paths from s to T. 
Then n = mincer is the max-flow min-cut capacity of 
the multicast connection. 

B. Network code 

We follow [3] in describing the network code. For each node 
w G V, let the set of all incoming edges be denoted by r/(u). 
Then |r7(w)| = Si{v) is the in-degree of v. Similarly the set 
of all outgoing edges is defined by Toiv), and the out-degree 
of the node v is given by |ro(w)| = 5o{v)- For any e G £ and 
w e V, let head{e) = v, if v is such that e G Tj{v). Similarly, 
let tail{e) = v, if v is such that e e To{v). We will assume 
an ancestral ordering on £ of the acyclic graph Q. 

The network code can be defined by the local kernel 
matrices of size Si{v) x 6o{v) for each node v <E V with 
entries from Fg. The global encoding kernels for each edge 
can be recursively calculated from these local kernels. 



The network transfer matrix, which governs the input-output 
relationship in the network, is defined as given in [3]. Towards 
this end, the matrices A,K,and _B^(for every sink T G T are 
defined as follows: 

The entries of the n x 151 matrix A are defined as 



A, 



where G F, 







if ej e Tois) 
otherwise 



q is the local encoding kernel coefficient at 



the source coupling input i with edge Cj G ro(s). 
The entries of the \£\ x \£\ matrix K are defined as 







if head{ei) — tail{ej) 
otherwise 



where the set of f3ij G is the local encoding kernel coeffi- 
cient between and ej, at the node v — head{ei) — tail{ej). 
For every sink T ^ T, the entries of the \E\ 'x n matrix 
are defined as 



eF„ 







if Gr,(r) 

otherwise 



where all eg 

For instantaneous networks, we have 

F {I-K)-^ 

where / is the \£\ x \£\ identity matrix. Now we have the 
following: 

Definition 5 ( [3]): The network transfer matrix, Mt, cor- 
responding to a sink node T G T is a full rank n x n matrix 
defined as Mt := AFB'^ = AFt- 

Definition |5] implies that if a; G F^' is the input to the 
instantaneous network at any particular instant, then at any 
particular sink T G T, we have the output, y G F^\ at the 
same instant, to be y = xMt- 

C. Convolutional codes for networks 

Assuming that a n-dimensional linear network code mul- 
ticast has been implemented in the network, we define the 
following terms- 

Definition 6: An input convolutional code, Cg is a convolu- 
tional code of rate k/n{k < n) with a input generator matrix 
Gi{z) implemented at the source of the network. 

Definition 7: The output convolutional code Ct, corre- 
sponding to a sink node T G T is the k/n{k < n) 
convolutional code generated by the output generator matrix 
Go,t{z) which is given as Go,t{z) = Gi{z)Mt, with Mt 
being the full rank network transfer matrix corresponding to 
a n-dimensional network code. 

Example 1: Consider the 4C2 combination F3 network 
as shown in Fig. [T] For this network, let the input con- 
volutional code over F3[z] be generated by G/(z) = 
[1 + 1 + z + z^] . The network transfer matrices at each 
sink and their corresponding output convolutional matrices are 
calculated and tabulated in Table |T] 

Thus, as can be seen from the Example [T] the source 
node implements a convolutional code, and maps the encoded 




Fig. 1. 4C2 combination network over a ternary field. Tlie global kernels of 
the edges coming from the source are indicated. All the intermediate nodes 
have local kernels unity. 

TABLE I 

4C2 - F3 NETWORK FOR THE INPUT CONVOLUTIONAL CODE 
Gl{z) =[l + Z^ 1 + 2 + 22] . 



Sink 


Network transfer 
matrix 


Output convolutional code 


Ti 


Mt, = 






Go,Ti(z) = [1 + 2^ 1+2 + 2^] 


T2 


Mt^ = 


( ^ ^ \ 




Go,T^{z) = [l + z^ 2 + 2 + 2^2] 


T'i 


Mt^ = 


M 1 1 

\Q 2 ^ 




Go,tA^) = [1 + 2^ 22] 




Mt^ = 


( Q \\ 
\^ \ 1 ^ 




Go,tJz) = [1 + z + 2 + 2 + 22^] 


n 


Mt, = 


^ 1 N 




Go,TrJz) = [1 + 2 + 22 2z] 


Te 


Mt, = 


I 1 2 J 




Go,Te{z) = [2 + 2 + 222 2z] 



symbols into its outgoing symbols. The network maps these 
symbols from the source to symbols at the receivers. Each of 
the receivers hence sees a different convolutional code which 
might have different distance properties and different degrees 
(S). Fig. ID illustrates the entire system for a particular sink. 

D. Network-errors 

An error pattern p, as stated previously, is a subset of £ 
which indicates the edges of the network in error An error 
vector It) is a 1 X |£ I vector which indicates the error occurred 
at each edge. An error vector is said to match an error pattern 
(i.e w G p) if all non-zero components of w occur only on the 
edges in p. An error pattern set $ is a collection of subsets 
of £, each of which is an error pattern. Therefore we have the 
formulation as follows. 

Let a; G F" be the input to the network at any particular time 

\£\ 

instant, and let w d Fq he the error vector corresponding 
to the network-errors that occurred in the same particular 
instant. Then, the output vector, y G F^\ at that instant at 
any particular sink T can be expressed as 

y — xMt + wFt 



Convolutional encoder of the 
output convolutional code at sink T 
G^^(z) = G/z)AFB^ 



Source 



Input 
convolutional 
code 
encoder 
G,(z) 



Local kernels 
at source 
(Matrix A) 



Network 
(Matrix F) 



Matrix 
B 



SinkT 



JConvolutional 
code 
decoder 



5,(T) 



r 



Fig. 2. A network with a input convolutional code and a network code 



IV. Convolutional Codes for network-error 
Correction - Code Construction and Capability 

A. Bounded distance decoding of convolutional codes 

In this section, we briefly discuss and give some results 
regarding the bounded distance decoding of convolutional 
codes, reinterpreting results from [18] for our context. 

For the convolutional encoder with c encoded output sym- 
bols and b input (information) symbols, starting at some state 
in the trellis, we shall denote every such c output symbol 
durations as a segment of the trellis of the convolutional code. 
Each segment can be identified by an integer, which is zero 
at the start of transmission and incremented by 1 for every c 
output symbols henceforth. 

Let C be a rate b/c convolutional code with a generator 
matrix G{z). Then corresponding to the information sequences 
Uo,Ui, ..{ui G F^) and the code sequence Vo,Vi, ...{vi G 
F^), we can associate an encoder state sequence cro,cri,.., 
where at indicates the content of the delay elements in the 
encoder at a time t.We define the set of j output symbols 
as i'[oj) ■— [vo,Vi, ...,Vj-i] . We define the set Sdf^^^ 
consisting of all possible truncated code sequences ■U[o.j) V 
j of weight less than dfree{C) that start in the zero state as 
follows: 



which thereby can be considered as a code property because 
of the fact that Sdf^^^ is invariant among minimal encoders. 
Then, we have the following proposition: 

Proposition 1: The minimum Hamming weight trellis de- 
coding algorithm can correct all error sequences which have 
the property that the Hamming weight of the error sequence in 
any consecutive Td^^^ (C) segments is utmost '^J'-"^'-) — ij . 

Proof: Without loss of generality, let at be a correct state 
(according to the transmitted sequence) at some depth t in 
the path traced by the received sequence on the trellis of the 
code and let us assume that all the errors before t have been 
corrected. 

Now consider the window from t to t+Td^^^^ (C), consisting 
of Td.^^{C) segments. In this window, the Hamming weight 



= (C)- 



— J . However, by the 



of the error sequence is utmost 

definition of Tdf^^^{C), the distance between the correct path 
and every other path of length Td^^^^ (C) starting from the state 
at is at least dfree{C)- Therefore, in this window, the error 
sequence can be corrected. 

Now using at+Ta (C) depth t + Tdf^^^{C) as our new 
correct starting state, we can repeat the same argument thus 
proving that the entire error sequence is correctable. ■ 



Sd 



I '^H ('^[Oj)) < dfree{C),ao = 0, V j > O} B. Construction 



where wh indicates the Hamming weight over F^. Clearly 
the definition of Sdj,^^^ excludes the possibility of a zero state 
in between (in the event of which wjf (■U[o.j)) ^ djree{C)), 
i.e, cTj ^ for any t such that Q < t < j. We have that the 
set Sdf^^^ is invariant among the set of minimal convolutional 
encoders. We now define 



T, 



-AC) - 



max 



J + 1 



For the given network with a single source that has to 
multicast information to a set of sinks, n being min-cut of 
the multicast connections, a ?i-dimensional network code in 
place over a sufficiently large field Fq (for which we provide 
a bound in Subsection lIV-Eb of characteristic p, we provide 
a construction for a convolutional code for correcting errors 
with patterns in a given error pattern set. This is the main 
contribution of this work. The construction is as follows. 



1) Let Mt — AFt be the n x n network transfer matrix 
from the source to any particular sink T T. Let $ 
be the error pattern set given. Then we compute the 
following sets. 

2) Let the set of all error vectors having their error pattern 
in $ be 

3) Let 

W-r := {wFt I w e W$} 

be computed for each sink T. This is nothing but the set 
of n-length resultant vectors at the sink T due to errors 
in the given error patterns p G $. 

4) Let 

Ws := U {w^M^^ I e Wt} 
TeT 

be computed. This is the set of all n length input vectors 
to the network that would result in the set of output 
vectors given by Wt at sink T, for each sink T. 

5) Given a vector y G F™ (for some positive integer m), 
let WH{y) denote the Hamming weight of y, i.e., the 
number of non-zero elements of y. Let 

ts = max wh{ws). (2) 

6) Choose an input convolutional code Cg with free distance 
at least 2ts + 1. 

C. Decoding 

Let Gi{z) be the k x n generator matrix of the input 
convolutional code, Cg, obtained from the given construction. 
Let Go,t{z) = Gi{z)Mt be the generator matrix of the 
output convolutional code, Ct, at sink T <E T, with Mt being 
its network transfer matrix. Each sink can choose between two 
decoding methods based on the free distance (d free{CT)) and 
Tdj^^^ (Ct) of its output convolutional code as follows: 

Case-A: This case is applicable in the event of both of the 
following two conditions are satisfied. 

dfree{CT) > 2 I max wh{w^) 1+1 (3) 

and 

Td,_(C,) >T<j,_(Ct). (4) 

In this case, the sink T performs minimum distance decoding 
directly on the trellis of the output convolutional code, Ct- 

Case-B: This case is applicable if either of the following 
two conditions are satisfied. 

d/ree(CT) < 2 ( max Wh{Wt)]+'^ 

or 



This method involves additional processing at the sink, i.e, 
matrix multiplication. We have the following formulation at 
the sink T. Let 

Wliz) V'^iz) ... V'^{z)] = [Vl{z) V2{z) ... Vn{z)] + 

[Wi{z) W2{z) ... Wn{z)\ 

represent the output sequences at sink T, where 

[vi{z) V2{z) ... Vn{z)] ^ u{z)Go,t{z) ^ u{z)Gj{z)Mt 
u{z) being the k length vector of input sequences, and 

[wi{z) W2{z) ... Wn{z)] 

represent the corresponding error sequences. Now, the output 
sequences are multiplied with the inverse of the network 
transfer matrix Mt, so that decoding can be done on the trellis 
of the input convolutional code. Hence, we have 

v'^iz) ... <(z)] = K(z) v'.iz) ... <(z)]M^i 

^ u{z)Gl{z) + [Wi{z) W2{z) ... Wniz)]M^^ 

= u{z)Gi{z) + [w[{z) ... w'^{z)] 

where w'{z) — [w'i{z) W2{z) ... ■wl-^{z)] now indicate the 
set of modified error sequences that are to be corrected. Then 
the sink T decodes to the minimum distance path on the trellis 
of the input convolutional code. 

D. Error correcting capability 

In this subsection we prove a main result of the paper 
given by Theorem [3] which characterizes the error correct- 
ing capability of the code obtained via the construction of 
Subsection IIV-BI Before proving the following theorem, we 
recall the following observation that in every network use, 
n encoded symbols which is equal to the number of symbols 
corresponding to one segment of the trellis, are to be multicast 
to the sinks. 

Theorem 3: The code Cg resulting from the construction 
of Subsection IIV-BI can correct all network-errors that have 
their pattern as some p G $ as long as any two consecutive 
network-errors are separated by Tdf^^^{Cs) network uses. 

Proof: In the event of Case-A of the decoding, the given 
conditions ((O and (|4]i) together with Proposition [T] prove the 
given claim that errors with their error pattern in $ will be 
corrected as long as no two consecutive error events occur 
within Tdf^^^{Cs) network uses. 

In fact, condition O) implies that network-errors with pat- 
tern in $ will be corrected at sink T, as long as consecutive 
error events are separated by T(;^^^^(Ct). 

Now we consider Case B of the decoding. Suppose that the 
set of error sequences in the formulation given, w'{z), is due 
to network-errors that have their pattern as some p G such 
that any two consecutive such network-errors are separated by 
at least Tdf,.^^{Cs) network uses. 

Then, from we have that the maximum Hamming 
weight of any error event embedded in w'{z) would be utmost 
ts, and any two consecutive error events would be separated by 
Tdf^^^{Cs) segments of the trellis of the code Cg. Because of 



the free distance of the code chosen and along with Proposition 
[T] we have that such errors will get corrected when decoding 
on the trellis of the input convolutional code. ■ 

E. Bounds on the field size and Tdf^^^{Cs) 

1) Bound on field size: The following theorem gives a 
sufficient field size for the required {n,k) convolutional code 
to be constructed with the required free distance condition 
(> 2is + 1). 

Theorem 4: The code Cg can be constructed and used to 
multicast k symbols to the set of sinks T along with the 
required error correction in the given instantaneous network 
with min-cut n (n > k), if the field size q is such that 

f 

n\q—l and q>max<\T\, - + 

n — k 

Proof: The condition that 

q>\r\ 

is from the known sufficient condition [4] for the existence of 
a linear multicast network code. 

For the other conditions, we first note that in the construction 
of Subsection IIV-BI tg < n. In the worst case that tg ~ n, 
we need dfree{Cs) > 2n + 1. We have from the generalized 
Singleton bound: 

dfreeiCs) < {u ^ k)([S/k\ +l)+5+l. 

In order that dfree{Cs) be at least 2n + 1, we let S = 2k, in 
which case the R.H.S of the inequality becomes 

{n-k){[2k/k\ +l) + 2fc+l 
= 2n + {n - k) + 1 > 2n + 1 

Thus, with S ^ 2k, from ^ we have that (n, k,5 ^ 2k) MDS 
convolutional code can be constructed based on [17] if 

2n2 „ 



n\q — 1 and q > 

11 — K 

Such an MDS convolutional code the requirements in the 
construction (dfree{Cs) > 2n + 1), and hence the theorem 
is proved. ■ 

2) Bound on Tdf^^^{Cs): Towards obtaining a bound on 
Tdf^^^{Cs), we first prove the following lemma. 

Lemma 1: Let C be a rate b/c convolutional code with 
degree 5 and Sdf^^^ be defined as in Subsection IIV-AI for a 
minimal encoder (a controller canonical form realization [14] 
of a minimal basic generator matrix, Gmb{z), of C). Then any 



cannot have & zeros in 5 consecutive segments, i.e, at least 
one of Vi,Vi, Vi^g^i is non zero V < i < j — (5. 

Proof: Let the ordered Forney indices (row degrees 
of Gmf,(z)) be 1^1,1^2, ■■■,i^b — i^max, and therefore S be- 
ing the sum of these indices. Then a systematic generator 
matrix(G'sys(z)) for C that is equivalent to Gmb{z) is of the 
form 

Gsysiz) = T~'^{z)Gmb{z) 



where T{z) is a full rank 6x6 submatrix of Gmb{z) with a 
delay-free determinant. We have the following observation. 

Observation 1: The degree of det {T{z)) is clearly utmost 
5. Also, we have the element ti,j{z) of T^^(z) as 



U^jiz) 



Cof actor {T{z)jj) 



det (r(z)) 

where Cof actor {T{z)j^i) G F^[z] is the cofactor of the 
element of T{z). The degree of Cof actor {T{z)j,i) 
is utmost 5 — Vj < 5 — vi. 

Let aij{z) G Fq(z) represent the {i,jY'^ element of 
Gsys{z), where 



J {z) 



J2iiMz)9k,jiz) 

k=l 



_ J2k=i Cofactor{T{z)k,i)gk,jiz) 
det{T{z)) 

gk.j{z) being (fc,j)*'' element of Gmb{z)- Therefore, the 
element aij{z) can be expressed as 

P^,jiz) 



a-,j{z) 



det (T(z)) 

where the degree of Pi,j{z) G ^q[z] is utmost S + Vr, 
Now if we divide Pi,j{z) by det {T{z)), we have 



aij{z) = qt,j{z) + 



n^] jz) 
det {T{z)) 



(5) 



where the degree of qij{z) e ^q[z] is utmost 



vi, and 



the degree of ri_j{z) is utmost 5 — 1. Because every element 
of Gsys{z) can be reduced to the form in (|5]l, we can have 
a realization of Csys{z) with utmost 6 memory elements for 
each of the 6 inputs. Let this encoder realization be known as 
E. 

Now we shall prove the lemma by contradiction. Suppose 
there exists a codeword v{z) = [vo,Vi, ...,Vj^2,Vj-i,Vj, ...] 
exists such that t'fo.j) = [va,Vi, ...,Vj^i] G Sdf,.^^ and 
Vi,Vi, ...,Vi^s_i are all zero for some i such that < i < 

j-s. 

Let Ug{z) be the information sequence which when encoded 
into v{z) by the systematic encoder E. Because of the 
systematic property of E, we must have that Ui, Ui, Ui+s-i 
are also all zero. By Observation [T] i? is an encoder which 
has utmost 6 delay elements (for each input), and hence the 
state vector (Ti+s at time instant i + 6 becomes zero as a result 
of these 6 zero input vectors. Fig. |3] shows the scenario we 
consider 

Therefore the codeword v{z) can be written as a unique 
sum of two code words v{z) = v'{z) + v"{z), where 

i+5-l 

v'{z)^ ^ Vkz^ ^[v^,...,Vi^O,...,v^+s^i^Q,0,..] 



and 

v"{z) 



k=l 



E 

k=i+S 



Vkz^ = [0, 0, 0, Vi+s, Vi+s+i, ...,Vj, ...] 



Number of segments = j 



Time t = 




If the Hamming weight of the path 
in all the consecutive 5 segments is 
zero, then a = 0. 



Fig. 3. The trellis corresponding to a systematic encoder of C 



where G and the uniqueness of the decomposition holds 
with respect to the positions of the zeros indicated in the two 
code words v'{z) and v"{z). 

Let Umb{z) be the information sequence which is encoded 
into v{z) by a minimal realization E„ib of a minimal basic 
generator matrix Gmb{z) (a minimal encoder). Then we have 



Umb{z) 



'•mb 



where u'^^^{z) and u'^f^{z) are encoded by E^ab into v'{z) 
and v"{z) respectively. 

By the predictable degree property (PDP) [14] of minimal 
basic generator matrices, we have that for any polynomial code 
sequence v{z), 

deg {v{z)) = max {deg {umbA^)) + ^i} ■ 

l<l<b 

where Umb.i{z) G ^q[z] represents the information sequence 
corresponding to the Z*'* input, and deg indicates the degree 
of the polynomial. Therefore, by the PDP property, we have 
that deg {u'^^^{z)) < i, since deg{v'{z)) < i. 

Also, it is known that in the trellis of corresponding to 
a minimal realization of a minimal-basic generator matrix, 
there exists no non-trivial transition from the all-zero state 
to a non-zero state that produces a zero output. Therefore we 
have deg {u'^i^{z)) > i + 5, with equality being satisfied if 
Vi+s 7^ 0. Therefore, Umb{z) is of the form 



u„ib{z) 



,{z) 



,{z) 



k=l 



k=i+S 



It k 

%b,kZ 



U,nb{z) — [M„,f,_Q, .., M„jj 0, 0, ..] 

+ [O, .., 0, M^f„i+5: '"^ 



mb^i+5+1 1 ' 



i.e. 



'^mb 



(z) 



y^mbjO^ '^rnb.l 7 • • ■ ; '^rtib/i 



) '^vab,i+5— 1 ; '^rabA+5 ; ■ ■] 



where UmbA = Umb,i+i = ■■■ = u„iba+6-i = e F^. 

With the minimal encoder E„ib, which has utmost i/f, 
memory elements, these 6 consecutive zeros of would result 
in the state vector (Tmb.t becoming zero at time instant 
i + i'b<'i + S<j, i.e, (Tmb.i+u^ = 0. But the definition of Sdf^^^ 
excludes such a condition, which means that ■U[o.j) ^ ^dj,^^^, 
contradicting our original assumption. Thus we have proved 
our claim. 

■ 

We shall now prove the following bound on Tc;^^„(C). 
Proposition 2: Let C be a (c, b, 6) convolutional code. Then 



{C)-l)S+l. 



(6) 



Proof: Let ^[oj) £ Sdf^^^ be some truncated codeword. 
Then we have wh (■^[o.j)) dfree{C) — 1. By Lemma [T] 
we have that in any consecutive 5 segments, the Hamming 
weight of I'fo.j) is at least 1. With this observation, and by 
the definition of Tdj^^^{C), we have (|6]l, thus proving the 
proposition. ■ 
Thus, for a network error correcting MDS convolutional code 
Cs, we have the following bound on Tdf^^^Cs). 

Corollary 1: If the code Cg chosen in the construction of 
Subsection IIV-BI is an (n, k) MDS convolutional code, then 
we have the following bound on Tdf^^^{Cs)- 



Tdj^^ACs) < 6nk - 2fc2 + L 



(7) 



Proof: In the Construction of Subsection llV-BI if the code 
Cs selected is an MDS convolutional code, then we know 
from the proof of Theorem |4] that the degree being 5 — 2k 
satisfies the required error correcting capability. Moreover, for 
an {n, k, 6) MDS convolutional code, we have 

dfreeiC) = (n - k){l6/k\ + I) + S + 1 

Therefore, substituting this value for dfree{Cs) with 5 — 2k 
in ^ of Proposition m we have (|7]i. ■ 

V. Illustrative Examples 
A. Code construction for the butterfly network 

The most common example quoted in network coding 
literature, the butterfly network, is shown in Fig. |4] Let us 
assume the ancestral ordering as given in the figure. Every 
edge is assumed to have unit capacity. It is noted that the 
network code in place for the butterfly network as shown is 
a generic network code for all field sizes. We seek to design 
a convolutional code for this network which will correct all 
single edge errors. 

Example 2 (Butterfly network under a binary field): The 
network transfer matrix for sink Ti is the full rank 2x2 
matrix 

" 1 1 
1 



= AFt 



where 



A^ 



100000000 
010000000 




B applies and decoding is done on the trellis of the input 
convolutional code after processing. 

Similarly, the convolutional code thus seen by the sink node 
T2 has the generator matrix 

,21 



Go,tM = Gi{z)Mt,{z) 



\z 1 + z + 



Fig. 4. Butterfly network 



1 

1 1 



Ft, = 

Similarly, for sink T2, A/jv 



1 
11110 



1 
1 1 



j4Ft, where 



Fto — 



1 1 
1 



1110 




1 
1 



For single edge errors, the error pattern set is 
* = {{ea:i = l,2,...,9}. 
Then the set of 9 length error vectors over F2, >V$ = 
{(1, 0,0,. .., 0),(0,1,0,...,0),..., (0,0,0,. ..,0,1), (0,0, 0,...,0)} 
For both sinks Ti and T2, we have 

Wt = Wt, = Wt, = {(0, 0), (0, 1), (1, 0), (1, 1)} 
Since A/^^^ = Mt, and = Mt2, we have 



This is a free distance 4 code, which is again lesser than 
2 (max.u,y^gvVr2 '^h{wt2)) +1 = 5. Hence, for this sink 
too. Case B applies and decoding is done on the input trellis. 

Example 3 (Butterfly network under a ternary field): We 
now present another example to illustrate the case when 
the field size and the choice of the input convolutional 
code affects the error correction capabilities of the output 
convolutional codes at the sinks. Let us assume the butterfly 
network with the network code being the same as the previous 
case, but over F3. The network transfer matrices in this case 
are the same as before, but the symbols are from F3. We seek 
to correct single edge errors in this case too. Thus the error 
pattern set is the same as the previous case. Now we have 
the set of 9 length error vectors over F3 

= {(1,0,0,.. .,0), (0,1,0,.. .,0),..., (0,0,0,.. .,0,1), 
(2,0,0,...,0),(0,2,0,...,0),...,(0,0,0,...,0,2),(0,0,0,...,0)} 

The other sets are as follows. Wt = = Wjb = 

{(0,0),(0,1),(1,0),(1,1),(0,2),(2,0),(2,2)} 



With AL 



Ti 



1 2 
1 



and Mrp ^ 

-1 2 



we have 



Ws= \J {w^Mt I e Wt} 



= {(0,0), (0, 1), (1,0), (1,2), (0,2), (2, 1), (2,0)} 
Thus, again we have, 

ts 



TeT 



max wh{ws) 



max wuiw^) — 2. 



W. = {(0,0), (0,1), (1,0), (1,1)} 



-1 = 



Hence a 

f 1 = 5 



Now we have ts = max^i^evVs '"^//(''^s) = 2 
convolutional code with free distance at least 2ts 
is required to correct these errors. With the min-cut n being 
2, let fc = 1. Let this input convolutional code Cs be generated 
by the generator matrix Gi{z) = \l + z"^ 1 + z + z^] . 

This code is a degree 2 convolutional code with free distance 
5, and Tdj^^^ {Cs) — 6. Hence, by Theorem |3] this code 
will correct all single edge errors under the condition that 
consecutive single edge errors occur within 6 network uses. 
Now the sinks must select between Case A and Case B for 
decoding, based upon their output convolutional codes. 

The output convolutional code that is 'seen' by the sink Ti 
has a generator matrix 

Go,t,{z) = Gi{z)Mt, = [1 + z' z]. 



Hence a convolutional code with free distance at least 2ts 
5 is required to correct all single errors. 

We compare the error correction capability of the output 
convolutional code at each sink for two input convolutional 
codes, Cs and C's generated by the matrices 
Gi{z) = [l + z^ i4.^^^2j 



and 



G;(z) = [1 + z2 1 + z + 2z2 



respectively, each over F3[z]. Both of these codes are degree 
2 convolutional codes and have free distances dfree{Cs) = 

dfreeiC's) = 5, with Td^^JCs) = Td.^^M) = 6- 

First, we discuss the case where the input convolutional 
code is Cs- The sink Ti thus sees the code generated by 

Go,Ti(z) = Gi{z)Mt, = [1 + 2 + Z + 2z^] 



This code seen by sink Ti has a free distance of only 3, which which has a free distance of 5, with Tdj,^^^{CTi) = 6 
is lesser than 2 (max^^^ gw,,^ WHiwTi)) +1 = 5. Hence case 



Tdf^^^ (Cs). Thus decoding is done on the output trellis at sink 



Ti to correct all single edge errors that as long as they are 
separated by 6 network uses. Sink T2 sees the code generated 
by 

Go,T2 (2) = [2 + z + 2z2 l + z + z^] 

which has dfree = 6, and ^^^.^^^(CTa) = 6 = Td^^^^JCs). 

Therefore, sink T2 can also decode on the output trellis after 
multiplication by the corresponding processing matrix to get 
the required error correction in every 6 network uses. Upon 
carrying out a similar analysis with the input convolutional 
code being C'^, we give the following tables for comparison. 



TABLE II 

Butterfly NETWORK WITH Cs[d/ree(Cs) = 5,Ttj^.^^^(Ca) = 6] 



Sink 


Output convolutional 
code [Go,Ti(z)] 


dfree{CTi ), 


Decoding on 


Ti 


[l + z^ 2 + z + 2z^] 


5,6 


Output trellis 


T2 


[2 + z + 2z-^ 1 + z + z'-'] 


6,6 


Output trellis 



TABLE III 

Butterfly NETWORK WITH C5[d/ree(Cs) = 5,Td^.^_,_,(C^) = 6] 



Sink 


Output convolutional 
code [Go,Ti{z)] 


dfreeiCTi), 


Decoding on 


Ti 


[1 + z^ 2 + z] 


4,3 


Input trellis 


T2 


[2 + z 1 + z + 2z^] 


5,5 


Output trellis 



With the input convolutional code being Cg, conditions (O and 
dU are satisfied at both sinks. Hence additional processing can 
be avoided at both sinks and they can decode on the output 
trellis directly, and get single edge error correction under the 
constraint that consecutive single edge errors are separated by 
at least 6 network uses. 

However with C^, one of the sinks Ti does not have 
sufficient free distance at its output convolutional code, and 
hence has to process the incoming symbols using A/jT^ and 
thereby decode on the trellis of the input convolutional code. 

Thus it can be seen that using a larger field size and 
choosing the input convolutional code appropriately can give 
more desirable properties to the output convolutional codes at 
the sinks. 

B. Code construction for the 4C2 network 

Example 4 (iC2 combination network under V'i): Let us 
consider the combination network in Fig. [T| The network 
transfer matrices for the 6 sinks are as in Table |T] We 
seek to design a convolutional code that will correct all 
network-errors whose error vectors have Hamming weight 
utmost 2 (i.e single and double edge errors). 

The eiTor pattern set is thus 

$ = {{e„ej}:i,j = l,2,...,15,16andf7^j} 

The set >V$ is the set of all 16 length vectors with Hamming 
weight utmost 2. We have 

Wt^ ^Wt2 = ... = Wt6 

and 

Ws= [j {w^Mj^^ I e Wt} = ¥l 



For every sink Ti, we have 

max wh(wt ) = max wh(w^) ~t^~2 

Therefore the input convolutional code needs to have free 
distance at least 5. 

As in Example [T] let the input convolutional code, Cs, over 
F3 [z] be generated by the matrix 

Gi{z) = [l + z^ l + z + z^]. 

This code has free distance = 5, and Tdj^^^{Cs) = 6. 

Each sink decodes on either the input or the output trel- 
lis depending upon whether dfree{CTi) > 2ts + 1, and if 
Tdfr-^^iCs) > Tdf^^^{CTi), and hence can correct all network- 
eiTors with with their pattern in $ as long as consecutive errors 
are separated by 6 network uses. The output convolutional 
codes at the sinks, their free distances and their Tdf^^^{CTi) 
are shown in in Table |IV] 



TABLE IV 

4C2 NETWORK WITH G/ (2) = [l + 1 + Z + z'^] 



Sink 


Output 
convolutional code 


dfreeiCTi), 
Tdf^„{CT^) 


Decoding on 


Ti 


[1 + 2^ I + Z + Z'-"] 


5,6 


Output trellis 


T2 


[1 + z^ 2 + z + 2z^] 


5,6 


Output trellis 


T-i 


[1 + z'^ 2z] 


3,4 


Input trellis 


Ti 


[1 + z + z'-' 2 + z + 2z^\ 


6,6 


Output trellis 


Ts 


[1 + Z + z'' 2z\ 


4,5 


Input trellis 




[2 + z + 2z^ 2z\ 


4,5 


Input trellis 



VI. Comparison with block network error 

CORRECTION CODES 

The approach of [7] can also be used to obtain network 
error correcting codes that correct t edge errors once in every 
J network uses (for some positive integer J). A time-expanded 
graph would then be used, i.e, with the network nodes (those 
except the source and sinks) and edges replicated for each 
additional time instant. 

Suppose the network has been replicated J times. Then the 
algorithm in [7] can be employed to obtain a t-error correcting 
BNECC for the time-expanded network, which equivalently 
for the original network gives a network error correcting code 
that corrects t errors once in every J network uses. It is noted 
that the sufficient field size q required by the technique of [7] 
to construct a i-error correcting BNECC for the time-expanded 
graph (T being the set of all sinks) is such that 




Our approach demands a field size according to Theorem ID 
which is independent of the number of edges in the network. 
Although the error correcting capability might not be com- 
parable to that offered by the BNECC, the reduction in field 
size is a considerable advantage in terms of the computation 
to be performed at each coding node of the network. Also, the 
use of convolutional codes permits decoding using the Viterbi 
decoder, which is readily available. 



For example, one could design network error correcting 
codes according to [7] for the butterfly network by using the 
twice replicated butterfly network as shown in Fig. |5] The 
time-expanded network has min-cut 4, and thus the technique 
in [7] can be used to obtain BNECCs, which correct single or 
double edge errors in the butterfly network once in 2 network 
uses. In either case, the sufficient field size q is such that 
q > 306, although by trial and error a code could be found 
over a smaller field size. On the other hand, the convolutional 
code that we used here in our paper for the butterfly network 
is over the binary and ternary fields. 



^1 




Fig. 5. A twice replicated butterfly network. Tlie edges are marked witli a 
time index as to denote tlie time-expanded nature of the network. 

VII. Discussion 

In the construction of Subsection IIV-BI the maximum 
Hamming weight tg of the vectors in the set Ws, is such 
that tg < n. Clearly the actual value of tg is governed by 
the network code and hence the network code influences the 
choice of the network-error correcting convolutional code. 
Therefore the network code designed should be such that tg is 
minimal, so that the free distance demanded of the network- 
error correcting convolutional code in the construction of 
Subsection IIV-BI is minimal. 

Also, for a particular error pattern set, the decoding proce- 
dure at the sinks (Case-A or Case-B of decoding as in Sub- 
section IIV-Q is influenced by the field size, the network code 
and the network-error correcting convolutional code chosen. 
The examples given in Section |V] illustrate the construction 
of Subsection IIV-BI and also compare the effects of change in 
field size and the convolutional code chosen to correct errors 
corresponding to a given fixed error pattern set. 

Acknowledgment 

This work was supported partly by the DRDO-IISc program 
on Advanced Research in Mathematical Engineering through 
a research grant to B. S. Rajan. 



References 

[1] R. Ahlswede, N. Cai, R. Li and R. Yeung, "Network Information 
Flow", IEEE Transactions on Information Theory, vol.46, no. 4, July 
2000, pp. 1204-1216. 

[2] N. Cai, R. Li and R. Yeung, "Linear Network Coding", IEEE Trans- 
actions on Information Theory, vol. 49, no. 2, Feb. 2003, pp. 371-381. 

[3] R. Koetter and M. Medard, "An Algebraic Approach to Network 
Coding", IEEE/ACM Transactions on Networking, vol. 11, no. 5, Oct. 

2003, pp. 782-795. 

[4] Ho, T. Medard, M. Koetter, R. Karger, D.R. Effros, M. Jun Shi 
and Leong, B., "A Random Linear Network Coding Approach to 
Multicast", IEEE Transactions on Information Theory, vol. 52, no. 10, 
Oct. 2006, pp. 4413-4430. 

[5] Raymond W. Yeung and Ning Cai, "Network error correction, part 1 
and part 2", Comm. in Inform, and Systems, vol. 6, 2006, pp. 19-36. 

[6] Zhen Zhang, "Linear network-error Correction Codes in Packet Net- 
works", IEEE Transactions on Information Theory, vol. 54, no. 1, Jan. 
2008, pp. 209-218. 

[7] Shenghao Yang and Yeung, R.W., "Refined Coding Bounds for network 
error Correction", IEEE Information Theory Workshop on Information 
Theory for Wireless Networks, July, 2007, pp. 1-5. 

[8] R. Koetter and F.R. Kschischang, "Coding for Errors and Erasures in 
Random Network Coding", IEEE Transactions on Information Theory, 
vol. 54, no. 8, Aug. 2008, pp.3579-3591. 

[9] D. Silva,F. R Kschischang, and R. Koetter, "A Rank-Metric Approach 
to Error Control in Random Network Coding", IEEE Transactions on 
Information Theory, vol. 54, no. 9, Sept. 2008, pp. 3951-3967. 
[10] E. Erez and M. Feder, "Convolutional network codes", IEEE Interna- 
tional Symposium on Information theory, Chicago, Illinois, USA, June 

2004, pp. 146. 

[11] N. Cai, R. Li, R. Yeung, Z. Zhang, "Network Coding Theory", 
Foundations and Trends in Communications and Information Theory, 
vol. 2, no.4-5, 2006. 

[12] S. R. Li and R. Yeung, "On Convolutional Network Coding", IEEE 
International Symposium on Information Theory, Seattle, Washington, 
USA, July 2006, pp. 1743-1747. 

[13] C. Fragouli, and E. Soljanin, "A connection between network coding 
and convolutional codes", IEEE International Conference on Commu- 
nications, Paris, France, 2004, vol. 2, pp. 661-666. 

[14] R. Johannesson and K.S Zigangirov, Fundamentals of Convolutional 
Coding, John Wiley, 1999. 

[15] G. D. Forney, "Bases of Rational Vector Spaces with applications to 
multivariable linear systems", SIAM J. Contr., vol. 13, no. 3, 1975, pp. 
493-520. 

[16] J. Rosenthal and R. Smaradanche, "Maximum Distance Separable 

Convolutional Codes", Appl. Algebra Engrg. Comm. Comput., vol. 10, 

no. 1, June 1999, pp. 15-32. 
[17] J. Rosenthal, H. Gluesing-Luerssen, and R. Smaradanche, "Construction 

of MDS convolutional codes", Appl. Algebra Engrg. Comm. Comput., 

vol. 47, no. 5, July 2001, pp. 2045-2049. 
[18] J.M.Jensen and I. S. Reed, "Bounded distance coset decoding of 

convolutional codes", lEE Proceedings, vol. 133, no. 5, pp. 1046-1060, 

1986. 



